Fast optimization for aircraft descent and approach trajectory 
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ABSTRACT 


We address problem of on-line scheduling of the aircraft de- 
scent and approach trajectory. We formulate a general mul- 
tiphase optimal control problem for optimization of the de- 
scent trajectory and review available methods of its solution. 
We develop a fast algorithm for solution of this problem using 
two key components: (i) fast inference of the dynamical and 
control variables of the descending trajectory from the low di- 
mensional flight profile data and (ii) efficient local search for 
the resulting reduced dimensionality non-linear optimization 
problem. We compare the performance of the proposed algo- 
rithm with numerical solution obtained using optimal control 
toolbox General Pseudospectral Optimal Control Software. 
We present results of the solution of the scheduling problem 
for aircraft descent using novel fast algorithm and discuss its 
future applications. 


1. INTRODUCTION 


In the future air traffic management system, the trajectory 
becomes the fundamental element of a new set of operating 
procedures collectively referred to as trajectory-based oper- 
ations (TBO) (Cate, 2013). The basis for TBO is that each 
aircraft’s expected flight profile and time (or airspeed) infor- 
mation will be specified by a four-dimensional (4D) trajec- 
tory (K. H. Shish et al., 2015; Young et al., 2016; K. Shish et 
al., 2016). 


One of the challenges in development of the future TBO is 
management of the airport congestion especially under con- 
vective weather conditions (M. Kamgarpour, W. Zhang, & 
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C.J. Tomlin, 2011). One of the key ingredients to the so- 
lution of this problem is fast online optimization of the air- 
craft descent trajectory. The difficulty in solving this problem 
stem from the fact that there are multiple phases (configura- 
tions) that have to be flown during descent while respecting 
the system dynamics and satisfying a large number of linear 
and nonlinear constraints. 


The standard approach procedure required by the ATC im- 
plies starting a continuous steady descent from 6,000 ft, or 
higher, following a steep descent to a set of cleared altitudes 
and joining the 3° glide-slope approach from below. The 
speed during final approach is based on the reference speed, 
V;ef, which is calculated on the basis of reference speed for 
Flaps 30 (for a specific type of the aircraft) and depends on 
the mass of the aircraft. 


In addition, the trajectory planning operation specifies a set 
of transitions to the given flaps and landing gear configura- 
tions with corresponding reference speed. The various phases 
flown by the aircraft during descent are controlled by the pi- 
lots and involve a set of nominal actions that are required to 
enable the designed vertical flight profile including e.g. cap- 
turing localizer and gliding slope. 


Overall, it can be seen that the flight planning problem render 
itself as a complex multiphase trajectory optimization prob- 
lem subject to dynamical, path, and control constraints (Betts 
& Cramer, 1995; Tomlin, Lygeros, & Sastry, 2000; Tomlin, 
Mitchell, Bayen, & Oishi, 2003; de Jong, 2014; de Jong et 
al., 2015, 2017). 


The resulting trajectory optimization problem can be ad- 
dressed using multiphase optimal control techniques. How- 
ever, the duration of the final approach is only several min- 
utes and changes in the flight plan may require fast online op- 
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timization. Current implementation of the optimization algo- 
rithm in an experimental Flight Management System (FMS) 
developed using General Pseudospectral Optimization Soft- 
ware (GPOPS) require more than 30 sec of optimization time 
from the top of descent to the runway (de Jong, 2014; de Jong 
et al., 2015, 2017). 


Furthermore, future implementations of the FMS demand 
continuous online estimations of the feasible time bounds for 
each phase of the flight and fast rescheduling of the flight plan 
in e.g. convective weather conditions. These demands call 
for analysis of alternative approaches to the online solution 
of multiphase trajectory optimization problem. 


Here we present a novel approach to the solution of this prob- 
lem and compare performance of our algorithm with the per- 
formance of a conventional technique. 


The paper is organized as follows. In the next section we 
provide formulation of the multiphase vertical trajectory op- 
timization problem. In Sec. 3 we analyze performance of the 
solution of this problem using General Pseudospectral Op- 
timization Software. In Sec. 4 we describe novel algorithm 
and its application to the optimization of the final approach 
to a runway in San Francisco airport. Finally, in Conclusions 
we summarize the obtained results and discuss future appli- 
cations of the algorithm. 


2. OPTIMIZATION OF THE AIRCRAFT DESCENDING 
TRAJECTORY 


2.1. Model Equations 


The system state in the vertical plane during the approach is 
defined as 


x ={V,7,2,h}, (1) 


where V is the speed, ¥ is the flight path angle, h is the alti- 
tude, and z is the distance to the runway. 


In addition, the aircraft control is represented by two virtual 
control inputs - thrust 7’ and angle of attack a 


u = {a,T}. (2) 


To simplify analysis without loss of generality we neglect 
mass change due to fuel burned and the wind. The resulting 
model takes the form 


mV =T cosa — D— mg siny, 


mV¥ = (Tsina +L) — mg cos 4, 3) 


ée =Vcos 7,he = V sin 7, (4) 


We note that actual control variable is pitch rate. However, 
due to time separation between slow and fast aircraft dynam- 
ics, the angle of attack a is considered to be virtual control 
input. 


The angles and forces used in the model are _ illus- 
trated (Miquel & Suboptimal, 2015) in Fig. 1 


0 Xx 


Figure 1. The aircraft forces and angles during vertical flight 
with nonzero climb rate. 


The lift and drag coefficients in the model are 
1 
p= 5ev?s (Cp, + Co,a+ Cp,,07+ 


Cds., be + Cp, Of + Ci Sig) ’ (5) 


1 
L= evs (Cr. + Cr,at+ C1,,, e+ 
itive be). 


where p is the air density, S is the net wing surface area, 65 
is the spoiler deflection, dy, is the flap deflection, di, is the 
gear setting. 


2.2. Problem formulation 


In general, the descent trajectory optimization requires the 
following problem to be solved (Betts & Cramer, 1995; Tom- 
lin et al., 2000, 2003; Becerra, 2010; de Jong, 2014; Pat- 
terson & Rao, 2015; de Jong et al., 2015, 2017). Find 
the control trajectories, u(t), t € [¢6?, w, state trajecto- 
ries c(t),t € ae | static parameters p\), and times 
©, on where 7 = 1,..., NV, is the number of phases, to min- 
imize the following performance index (Becerra, 2010): 


Np 
P< y (° [2 (e),p®, i) cs 
t=1 (6) 


subject to the dynamical constraints: 


#O() = f [eu Or, 4] 2 HD .EPL, 
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the path constraints 
aD <r (2), ue), 4] <a,  ® 
fort € e, i, the event constraints: 
e < e®) fag sue By eo) the | S el, (9) 


the phase linkage constraints: 


Wi< Va? ue (1 Oot ae 


: (10) 
aXe) up”, ae ule), pe), t0Ne) Saal <wU, 
where al) a) Ge ‘) y) and uf, he Mag, ) 
the bound constraints: 
(4) ee mG (t) < ul), te (eed 
py? < p® <p, (11) 
eee vase a 
(i) () & 72) 
ty a 
and the time constraints: 
2 >. (12) 


In each phase functions and variables in (6)-(11) have appro- 
priate dimensions that may change from phase to phase. 


2.3. Flight phases 


The following problem will be considered. Optimize the final 
approach trajectory using objective function (6) for the de- 
scending flight between CEPIN and the stabilized approach 
fix (500 ft) of the RWY 28R in SFO. The phases included 
into the analysis are shown in the Table 1. 


This simplified schedule was developed by taking into ac- 
count general requirements for the descent trajectory, see 
e.g. (Prats et al., 2014). Some additional limitations were 
also included. E.g. it was assumed that the LOC signal ex- 
tends 18 nm a way from and 4500 ft above the antenna site 
(see e.g. FAA-S-8081-9B, June 2001). It was also assumed 
that the glide slope capture is initiated 2 + 0.75 nm away 
from DUMBA in a level flight. In addition, acceleration and 
altitude increase during descent were excluded from this op- 
timization test. 


A number of optimization problems can be formulated within 
this framework. The critical parameters of interest during air- 
port congestions are earliest and latest time of arrival. A pa- 


rameter of common interest is fuel consumption. In this work 
we were primarily interested in scheduling transitions times 
between various descent phases and estimation of the earliest 
and latest transition time for each phase. 


Table 1. Phases included into the first test 


#]) Name | ALT GD) DST (am) | SPD (kn) 
1 | intial | 10000 29 230 
state 
alae | SOOO te | 292s we [a5 cy 
acu | < 10000 | ARCHI < 230 
; Seeeen 6000 < h Caen 2 sey 
ZILED | = 10000 | ZED ~ | 2230 
«| Pes | 3000 oS Gey 
P < 10000 Gree ~ | < 230 
5 | Loc | 4000 < h eae 2) t= 
capture < 4500 DUMBA < 230 
DUMBA 
@° | apis) eg | Sk ote ¥ 
= CEPIN = 
CEPIN 
7 | Flaps20 Bs z < a@< Ee Mi 
= AXMUL2 | © 
g | Gear | 1800 < h pea 165 < V 
down < 4500 AXMUL ~ < 195 
9 | GS 1800 < aaa 165 < V 
capture < 2000 AXMUL41 < 195 
AXMUL < 
iO: | ilaps2s+|| Pen |e & Axe || oY 
n(y£07) | muLs3 = 185 
DST x | AXMUL#2 | 150 <-V 
IL | Flaps30 | pn(yt6y) | < 2 <FIX | <170 
stabilized] 
12 | oom] 446 14 150 


We will now consider solution of this problem using freely 
available package Gauss Pseudospectral Optimization Soft- 
ware (GPOPS) (Rao et al., 2010) 


3. EXAMPLE OF NUMERICAL 
GPOPS 


SOLUTION USING 


The GPOPS package is easy to install and to use. It was 
shown to perform well for multiple aerospace applications in- 
cluding optimization of the descent aircraft trajectory at Na- 
tional Aerospace Laboratory (Netherlands). And we have 
chosen this package for initial evaluation for application to 
the scheduling problem. Because the convergence of the al- 
gorithm was quite slow we have initially chosen a subset of 
phases from the Table 1. 


Here we provide an example of optimization for descending 
trajectory with the following five phases 


1. flaps 15 gear up: V(2. = 215 kn; V1) = 95.3 kn; 


men 


2. flaps 15 gear down; Vie), = = 215 kn; V2) =953 kn; 


min 


ANNUAL CONFERENCE OF THE PROGNOSTICS AND HEALTH MANAGEMENT SOCIETY 2017 


3. flaps 20 gear down; V2), = 195 kn; v2) = 91.1 kn; 


4. flaps 25 gear down; Vie = 185 kn; vi = 87.5 kn; 
Note, that the method allows for many different types of ob- 
jective function (Rao et al., 2010). For example, the final 
altitude in each phase could be minimized or maximized us- 
ing functions similar to (14) and (15) and substituting time t¢ 
with altitude h. 


Alternatively, to minimize the total mechanical energy to fly 
along a given path the objective function can be chosen as 


ay 
J= ae V(t) - T(t) dt. (13) 


In the presented analysis the objective function of the opti- 
mization (the performance index in eq. (6)) was chosen to 
minimize 


— 4(1) (2) (3) (4) (5) 
J=t; +t +t +t + tp (14) 
or maximize 
= (1) (2) (3) (4) (5) 
J=—(ty +t +t +t +t ) (15) 


transition time between phases. 
The dynamical constraints in each phase t € e, i with 
1 = 1,...,5 (cf eq. (7)) were given by 


mV =T cosa — D— mg siny, 

mV¥ =(Tsina +L) — mg cos 4, 

te = V cos ¥, (16) 
he =Vsin Y, 


T =(T.—T)/rr, 


The bounds on the aircraft speed in each phase were given by 
Vion and V") 


min Usted above for each flaps configuration. 

The results of the minimization of transition times for four 
phases are shown in Fig. 2. It can be seen from the figure 
that transitions to the configurations gear down, flaps 20 and 
25 occur as soon as velocity of the aircraft approaches the 
corresponding limiting value. The flight path angle 7+ stays 
close to the lower bound corresponding to the fast descent. 
The convergence time was found to be very sensitive to the 
parameters of the problem and varies between 10 and 50 sec. 


The convergence time for maximization problem can vary be- 
tween 50 and 300 sec and the convergence is not robust. We 
note that the best reported performance of the optimized im- 
plementation of the multiphase pseudospectral algorithm (de 
Jong, 2014) in C++ using PSOPT package was 30 sec for the 
whole descent and final approach trajectory. 


We therefore conclude that although GPOSP package is po- 
tentially very useful for trajectory optimization at present it 
can only be used for off-line applications. To enable fast on- 


(a) (b) 
= 200 a -2 
ae one 
= = 
Ss 150 = -4 
100 -6 
(0) 100 200 te) 100 200 
t t 
3000 150 
(d)| (e) 
=. 2000 = 100 
= 1000 = 50 
te) oO 
(0) 100 200 te) 100 200 
t 


Ho 


at), deg 
ohRoNbOMORN 


Control: a(t) 


o so 100 150 200 
t 


Figure 2. Minimization of transition times for descending 
trajectory with 4 phases. The figures show the dependence 
of dynamical variables of the system on time: (a) Velocity; 
(b) y; (c) distance; (d) Altitude; (e) Thrust; and (f) Angle of 
attack. Different colors indicate different phases of the flight 
and correspond (left to right) to phases (i) to (iv). 


line optimization of the multiphase descending trajectory we 
propose a novel algorithm, which is considered in the next 
section. 


4, FAST MULTIPHASE OPTIMIZATION ALGORITHM 


Before we introduce the algorithm, let us provide some esti- 
mations for the final ILS approach to RWY 28R at SFO. The 
starting point for our analysis is the transition from fast to 
slow deceleration, which may normally happen at the flight 
level FL100 and speed Vo4g < 250 kn. 


According to the ”3:1 rule of descent” this transition point is 
located ~33 nm from the runway. At the point of stabilized 
approach the aircraft speed should be ~140 kn. So during 
descent the altitude and speed are reduced by approximately 
330 ft and 3.3 kn per each nautical mile. 


In addition, the aircraft has to capture the 3° gliding slope 
from below at the altitude approximately 1800 ft and distance 
~5.5 m from the runway, which translates into ~740 ftm ver- 
tical speed at the ground speed 140 kn. 


The flight plan normally should also accommodate a number 
of actions including: (i) localizer intercept; (ii) setting a se- 
quence of flaps configurations; (i11) deploying landing gear; 
(iv) capturing gliding slope; (v) initiating flare; (vi) chang- 
ing to touchdown phase. All this actions set additional con- 
straints on the flight profile. Note that there can also be mul- 
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ry, a, deg 
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Vdot [kn/s], gdot [rad/s] 
b 
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| 
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Figure 3. Aircraft dynamics reconstructed from the flight pro- 
file defined as h(x) and V(x): (a) a(t) and y(t); (b) V(t); (©) 
+(t) and V(t) (d) T(t) - thrust, D(t) - drag, and L(t) — W - 
lift minus weight; (e) flaps and landing gear configuration. 


tiple ATC corrections (constraints) to the normal approach 
procedure that have to be included into the flight profile. 


The large number of phases (aircraft configurations), dynami- 
cal constraints, and bounds on the control and dynamical vari- 
ables render trajectory optimization a complex multiphase 
optimization problem, see Table 1, cf. with 14 phases defined 
in (de Jong, 2014; de Jong et al., 2015, 2017). 


On the other hand, we see that a large number of phases 
must be accommodated on relatively short distance along the 
descending path. This fact can be used to substantially re- 
duce the dimensionality of the problem by approximating 
segments connecting neighboring transition points with low 
degree polynomials (in particular, straight lines). This as- 
sumption also agrees with the results of the optimization ob- 
tained GPOPS code, see Fig. 2 and cf. (de Jong, 2014; de 
Jong et al., 2015, 2017). 


Using polynomial approximation of the flight trajectories 
connecting phase transition points the dimension of the op- 
timization problem is reduced to 


(Nir — 1) x Par X Dip (17) 


Here N;, is the number of transition points between the 
phases of the flight, D;, is the dimension of each transition 
point (e.g. distance, altitude, speed), and P;, is the degree of 
the polynomial connecting two neighboring transition points. 
In this formulation the decision variables are the coefficients 


10000 


8000 


2000 


Figure 4. Bounds on the location of the events obtained us- 
ing Table | are shown by colored transparent parallelepipeds 
in 3D space (V - SPD, x - DST, h - ALT). The initial lo- 
cation of events is shown by open gray green squares (con- 
nected by green dashed line) located at the centers of the par- 
allelepipeds. The optimal solution is show n by the open blue 
circles connected by the solid blue line. 


of the polynomials and the locations of the phase transition 
points. 


Furthermore, using this approximation one can avoid dynam- 
ical constraints all together by reconstructing all the dynami- 
cal and control variables from the piece-wise polynomial ap- 
proximations of the flight profile. 


4.1. Inferring descending trajectory from vertical profile 


To infer the model dynamics from the flight profile we notice 
that eqs. (3) and (4) can be rewritten in the form 


ul = eosy (Tcos a — pSuCp — D—Wsiny) 
y= amas y (E sin a + pSuCz — W cos 7) (18) 
h’ =tany, 
using transformation 
dx = V cos(y)dt (19) 


and introducing distance x as a new independent variable 
along the flight path (Vinh, 1981). In eqs. (18) W = mg, 
u is the specific kinetic energy i and prime refers to deriva- 


tive with respect to x, e.g. h’ = dh/dz. 


Since V and h are known along a given flight profile, thrust 
T, angle of attack a, and flight path angle y can be found as 
functions of distance using eqs. (18) and then as functions 
of time using eq. (19). The results of calculations are shown 
in Fig. 3. In this example the flight profile corresponds to 
the spline approximation of the initial guess of the trajectory 
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corresponding to the location of transition points shown by 
the open gray circles in Fig. 4. 


It can be seen from the figure that all the dynamical and 
control variables can be successfully inferred from the given 
flight profile. The obtained results allow one to avoid colloca- 
tion methods (C.R.Hargraves & S.W.Paris, 1987) of dynami- 
cal trajectory optimization and obtain online solutions to the 
problems of primary interest for aircraft descent operations 
including minimization of additional drag and thrust during 
descent, and enforcing required time of arrival. 


To solve these problems we have to combine the algorithm, 
outlined above, with the optimization algorithm for the loca- 
tion of the transition points. 


4.2. Fast optimization algorithm 


Using piece-wise linear approximation to the flight profile 
and location of phase transitions points as decision variables, 
we reduce the problem of multiphase optimization of the de- 
scent trajectory to the following standard NLP: 


minimize f(x) 
x 


subjectto g;(z) <0, i=1,...,m, 
hj(z) =0, j=1,...,n, (20) 
A-x“<b, 
Aeq(2) < beg, 
lb<a<ub. 


Here vector of decision variables x has dimension N;, x Dt, 
number of phase transition pints times dimension of each 
point. 


The performance index (cost function) f(x) can have many 
different objectives. In the context of the descent trajec- 
tory optimization the most common objectives are minimiza- 
tion of the fuel use (~ minimum thrust) and required time 
of arrival (RTA) for maximum throughput of a given run- 
way/airport. 


In this work we do not consider the change of the phases or- 
der. The phases order will be fixed as shown in the Table 1. To 
enforce the phase order and no-climb, no-acceleration condi- 
tions we use inequality constraints in the form 


Begs tee, (21) 


where wx is the vector of decision variables {x1,..., Xer, hi, 


2+, er, Vi,---, Vir}, and the block-bidiagonal matrix Aeg is 
-1 1 0 0 
0 -1 1 0 
Aeg = ; 
0 0 -1 1 


nm 
Ss 
S 
= 
oO 
z= 
°o 
\ 
al 


7, @, deg 
V [knt], hdot [1 Oft/min] 
Ec 
(| 
Vdot [kn/s], gdot [rad/s] 
\ 


100 rT 
& s}/—v ]\ | 


5 0 15 
0 200 400 600 0 200 400 600 0 200 400 600 
time, sec time, sec time, sec 
200 30 400 
d e | - f 
ve OL nes | oya5/ 0) a | @ 
0 = 300 }\ 
Z > 20 2 % 
x —= ES eo X 
~ 100 T = 2 sot 
: —o O 15}) fas | 200 <—s 
Qa ——LW oD ———~ gear 7 \ vue 
) 50 : 2 ‘ 2 
7 a" = 10 ES 
0—— oo & 5 > SW i 
> ——<hh " 
-50 0 0 
0 200 400 600 0 200 400 600 3020 10 0 
time, sec time, sec DST, nmi 


Figure 5. Aircraft dynamics obtained as the result of opti- 
mization of the descending trajectory using local solver fmin- 


con (MATLAB): (a) a(t) and 7(¢); (b) V(t) and h(t); (c) ¥(t) 
and V(t) (d) T(t) - thrust, D(t) - drag, and L(t) — W - lift 
minus weight; (e) flaps and landing gear configuration; (f) 
speed and altitude as functions of the distance. 


By setting vector be, to zero we ensure that values of the dis- 
tance 2, altitude h, and speed V at a given step are no larger 
than the corresponding values at the preceding step. 


Bounds and the initial guess of the trajectory are shown in 
Fig. 4. It can be seen from the figure that bounds on the speed 
and altitude of the phase transition points can be both dis- 
connected and overlapping. It can also be inferred from the 
figure that the initial guess of the descending trajectory was 
obtained as a set of locations of the centers of the bounding 
boxes. 


In the first test we choose minimization of thrust as the objec- 
tive of the problem: 


f(z)= i- T° (a) dx. (22) 


0 


The solution of the NLP problem (20) - (22) was obtained us- 
ing a number of optimization solvers including local search 
with of fmincon (MATLAB), and global solvers such as 
genetic algorithm (MATLAB) and pattern search (nomad, 
OPTI) (Currie & Wilson, 2012). 


The best results were obtained using with local search fmin- 
con (MATLAB) as shown in Fig. 5. It can be seen from 
the figure that optimized values of the thrust are indeed very 
small, cf. Fig. 3. The corresponding variations of the flight 
path angle are also small, see Fig. 5(a). In addition, it can 
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be noticed from the figure (f) that the optimal vertical de- 
scent path is a smooth trajectory respecting the constraints. 
All these features are expected to be the main features of the 
descent trajectory with minimum thrust. 


We note that the vertical velocity is above 1000 ft/min for 
most time during descent. This is a quite high value and ad- 
ditional nonlinear constraints 


| V-sin(y) |< 700ft/min (23) 


may have to be imposed to keep vertical speed within prede- 
fined limits. 


Importantly, the convergence of the local search algorithm 
was consistently less than 5 sec. This result may suggest that 
there exist some strong, albeit hidden, convex properties of 
the objective function. 


For comparison, the convergence of the global search algo- 
rithms (genetic algorithm and pattern search) was slow. De- 
spite long convergence time the thrust found by the global 
search algorithms was much larger than the one described 
above. This fact also indicates that the cost function has 
strong convex property. 


Some insight into the properties of the cost function can be 


gained by rewriting it in the form (keeping terms ~ a?, y") 


i Powe SS ae (mv a Wy) 


pgs Uta)" ((¢ Co.) o (24) 
+ (S -Cp,,) a”) -Cp+we] , 


where V; are the mean value of speed in each phase, éV are 
the speed variations around the mean value, Cr = Cp, + 
Coy, Oft + Cp,, Sg: and Cp = Ch, + CLs Oft + CL,, St9- 


From the eq. (24) one can conjecture that the cost function is 
convex with respect to “hidden” optimization variables a and 
OV during the time interval At; corresponding to each phase. 
In addition, it is expected that there is a smooth continuous 
dependence of the thrust on the location of the boundaries of 
each phase. However, the vector of decision variables has 36 
components and the full analysis of the convex properties of 
the cost function is beyond the scope of this work and will be 
considered elsewhere. 


In this work we were primarily concerned with estimation of 
the earliest and latest transition times for each phase of the 
flight and accordingly the earliest and latest arrival times. To 
find these estimations the cost function was chosen as 


Ntr-1 


Cost = > At; 
1 


0.5 


a [deg] 


Ys 
Vdot [kn/s], gdot [rad/s] 
So 
a 
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Figure 6. The results of the minimization (dashed lines) and 
maximization (solid lines) of the required time of arrival. The 
notations are the some as in Fig. 5. 


to minimize arrival time and in the form 


Ntr-1 


Cost = — > At; 
i=1 


to maximize it. 


The solutions of the corresponding optimization problems is 
shown in Fig. 6. The obtained results show that the arrival 
times for the descent flight satisfying all the constraints are 
in the interval 484 + 600 sec. The time windows for the 
pilot actions during the descent are most clearly seen in the 
Fig. 6 (e) for the flaps and landing gear configurations. The 
corresponding changes in the flight profile can be observed 
in the Fig. 6 (f). It can be seen that the shortest time corre- 
sponds to a more steep descent with larger speed along the 
whole profile. These features are expected intuitively. They 
also confirm and quantify the intuitive idea that the aircraft 
speed during descent is the key scaling factor for the time 
windows allowed for the pilot actions. 


We note that these results were obtained without constraints 
on the vertical speed of the type (23). It is expected that con- 
straints will reduce the allowed time windows for pilot actions 
and for the arrival time. 


5. CONCLUSIONS 


We analyzed the problem of scheduling of the descent and ap- 
proach of commercial aircraft. The problem was formulated 
as a general multiphase optimal control problem. 
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We developed solution of this problem using Matlab package 
GPOPS. To evaluate the performance of this package for spe- 
cific applications to the scheduling problem we considered 
problem with 4 phases due to different flaps and landing gear 
configurations. It was shown that the package can solve com- 
plex multiphase problems in general form. However, its on- 
line applications are at present limited by slow convergence 
time. 


We proposed novel fast algorithm of scheduling and op- 
timization of the descent trajectory that reduces the origi- 
nal multiphase optimal control problem to the standard NLP 
problem of low dimension. It was shown that the dimension 
of the optimization problem can be reduced by 


e using algorithm of reconstruction of the full set of dy- 
namical and control variables along the descent path 
from a low-dimensional vertical and speed profile; 


e avoiding collocation methods of dynamical trajectory 
optimization; 


e choosing the location of the phase transition points as the 
decision variables; 


e approximating flight paths connecting transition points 
by low-dimensional polynomials. 


An additional advantage of the proposed algorithm is the abil- 
ity to include pilot actions during descent (such as captur- 
ing LOC and glide slop, changing flaps and gear configura- 
tion, establishing stabilized approach configuration, initializ- 
ing FLARE and touchdown etc.) directly into the flight plan. 


Preliminary testing of the algorithm shows promising results 
for the future online applications. In particular, the proposed 
approach paves the way to development of a general opti- 
mization algorithm that combines optimization of arrival time 
during descent of individual aircraft with optimization of run- 
way scheduling problem with multiple aircrafts, including sit- 
uation with convective weather conditions. 
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